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We present and implement a new scheme for independent control of both the final and the transient 
population utilizing the symmetry properties of the system. By proper pulse shaping, utilizing the 
invariance of the two-photon absorption to specific phase transformations of the pulse, different time 
evolutions of the transient population are photo-induced for a given (fixed) final state population. 
The model system is the Na atom. The work is conducted in the weak-field regime for which 
the transient two-photon excitation is described by second-order perturbation theory. One most 
attractive case is the extended family of third-order chirp pulses which control the population build- 
up duration independently of the final population. 



The general concept of coherent control is to steer a 
system towards a desired outcome by utilizing a manifold 
of interference pathways induced by a broad femtosec- 
ond pulse [TJ[31|3J|3]. Till now, this desired objective was 
mainly to control the transition probability to at least one 
of the given (fixed) final states. Such control was exper- 
imentally demonstrated in various atomic and molecular 

systems [sjiaHiiaiTniiiiJiiiiiiiiiiisiiisjiniiiHi 

IT]?! I201 |2"T] . However, we believe that complete control 
over the transient as well as final population is of ma- 
jor importance to any forthcoming application involving 
multiphoton processes and fast decoherence or relaxation 
rates. Despite the grate progress in coherent control of 
state-to-state transition probability, little work was done 
in controlling the transient population build-up. Most of 
the experimental advances in this area concentrated on 
experimental observation [22] and control [221 HI] of co- 
herent transients in a two-level atomic system, using one 
photon absorbtion process which does not allow control 
over the final population of the system, but only over its 
transient. Moreover, several important applications were 
also demonstrated, including femtosecond spectral elec- 
tric field reconstruction [55], high precision calibration 
of pulse-shaping setups [26 and most importantly, time 
resolved quantum state holography |27j . 

In this Letter we present for the first time experimental 
observation and full coherent control of transient state 
population excited via nonresonant two-photon transi- 
tion by shaped femtosecond pulses. The control scheme 
extends the frequency-domain picture [B] to describe the 
transient evolution of the population build-up. Then, 
by a proper pulse shaping, utilizing the invariance of 
the two-photon absorption to specific phase transforma- 
tions of the pulse [21], different transient evolutions of 
the population are photo-induced independently of a de- 
sired (fixed) final state population. The intuitive theoret- 
ical framework is also verified experimentally. Moreover, 
extended family of chirp pulses were found to be very 
robust in controlling separately the final and transient 
populations on a picosecond timescale using combination 
of frequency and time domain descriptions. The work is 



conducted in the weak-field regime for which the tran- 
sient two-photon absorption is described by second-order 
perturbation theory, which is the lowest order perturba- 
tion for final state phase control. The model system is 
the Na atom. In general, the essence of the present work 
is to provide a rational intuitive way to the full control 
over the state's population evolution as well as its final 
value on both, short and long timescales. 

We consider an atomic two-photon absorption process 
from an initial ground state |<7) to a final excited state 
|/), which is coupled via a manifold of states \v) that are 
far from resonance and have the proper symmetry. The 
light-matter interaction with a shaped temporal electric 
field e(t) is described by second-order time-dependent 
perturbation theory. Thus, the time-dependent ampli- 




FIG. 1: (a) Energy level diagram of the pump-prob excita- 
tion and detection scheme of Na. (b) Representation of the 
selective transient coherent control scheme including the base 
phase pattern and the anti-symmetric phase additions, (c) 
Pump pulse experimental spectra, (d) Outline of the experi- 
mental setup. 
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tude af(t) of state |/) at time t, is given by 
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with transition frequency f2: ^4o„_,. es is proportional to 

AW(uj fg ), while A^.^t) integrates all A^(u fg -S) 
with <5^0. The transient two-photon transition prob- 
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where /Ujj = (i| /i |j) is the dipole matrix element between 
a pair of states and u>ij = (Ei — Ej)/Ti is the correspond- 
ing transition frequency. Within the frequency-domain 
framework, the spectral field E(uj) = \E(uj)\ exp [z<I>(cj)] 
is given as the Fourier transform of e(t), with 
and <&(cl>) being the spectral amplitude and phase at fre- 

(2) 

quency uj. The transient amplitude A tr (t) can be written 
as 



In agrement with previous studies, by evaluating Eq. ^ 
for time t — > oo (i.e., after the pulse is over), one ob- 
tains A^ ] ear _ res {t oo) = A^_ res for finite ampli- 
tude of state |/), and the transition probability becomes 
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where w/ s and ^i 2 ^ are, respectively, the \g)-\f) tran- 
sition frequency and effective non-resonant two-photon 
coupling. The contribution to the transient amplitude 
can be separated into two channels of excitation: 
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The above frequency domain description allows intu- 
itively to identify the interfering pathways that coher- 
ently contribute to the transient two-photon transition 
amplitude (see Fig. 111. The time independent amplitude 

(2) — 

A OJ l_ res coherently interferes all the possible two-photon 
pathways from \g) to |/) (6 = 0), i.e., coherently in- 
tegrates over all their corresponding amplitudes. Each 
such pathway is composed of two absorbed photons of 
frequencies uj and ujf g —uj. As opposed to the time depen- 

(2) 

dent amplitude A ni ! ar _ res (t) which coherently interferes 
all the possible near-resonance two-photon pathways hav- 
ing a corresponding detuning (<5 ^ 0). It involves a non- 
resonant absorption of two photons with a two-photon 
transition frequency u>f g — 5 with a 1/(5 amplitude weight- 
ing and a time dependent phase term exp(— iSt). The on- 

(2) 

resonant pathways are excluded from A ni ; ar _ res (t) by the 
Cauchy's principal value operator p. The different am- 
plitudes are expressed using the parameterized amplitude 
A^ (O) interfering all the possible two-photon pathways 



Utilizing the symmetry property of A^ (uj f g ) that is 
invariant to any antisymmetric phase transformation, the 
two channels of excitation (i.e., the on-resonance time in- 

(2) 

dependent channel A or [_ res and the near-resonance tran- 

(2) 

'-isient channel A.^ ar _ res (t)) can be selectively controlled. 

IT 

Any phase ^* antisy m ^- ^ ^ na ^ ^ s antisymmetric around 
Wfg/2 keeps the value of A^ 2 \ujf g ) unchanged, while it 
generally changes the value of A^ 2 '(u}f g — 5) for 5 ^ 0. 

(2) 

Hence, it alters A y n ^, ar _ res (t) while keeping A y 01 [_ res in- 
variant. The effect of antisymmetric phase addition is ex- 

f 2) 

plained in details in [2T]. A i ^_ resl on the other hand, can 
be controlled from a zero till to its maximal value (TL- 
level) using a simple 7r-step spectral base phase ^^.^(a;). 

The total phase pattern = ^baseH + ^lntisym^ 
is presented schematically in Fig [T] 

The physical model system of the study is the sodium 
(Na) atom, with the 3s ground state as \g) and the 4s 
final state state as |/), (see Fig. [TJ. The transition fre- 
quency uifg = W4s.3s = 25740 cm -1 corresponds to two 
777-nm photons. The 3s-4s two-photon coupling is non- 
resonant via the p states. The Na atom interacts with 
phase-shaped linearly-polarized femtosecond pulses hav- 
ing a one-side blocked at 779 nm (12837 cm -1 ) Gaussian 
intensity spectrum centered around 790 nm (12658 cm -1 ) 
with 4 nm (66-cm -1 ) effective bandwidth (^200-fs TL 
duration), (see Fig. [IJ. 

The sodium vapor in a heated cell is irradiated with 
such laser pulses, after they undergo shaping in an op- 
tical setup incorporating a pixelated liquid-crystal spa- 
tial light phase modulator [28 . The effective spectral 
shaping resolution is Suj s haping=2-05 cm^ 1 (0.125 nm) 
per pixel. The peak intensity of the TL pulse is below 
10 9 W/cm 2 , corresponding to the weak-field regime. To 

(2) 

demonstrate the |/) state transient population P} r '(t), 
the 4s state is probed to the 7p state by weak unshaped 
linearly-polarized femtosecond pulses having a Gaussian 
intensity spectrum centered around 790 nm (12658 enr 1 ) 
with 12-nm (300-cm^ 1 ) bandwidth (~80-fs TL duration), 
at time delay t p . The experimental setup is presented 
schematically in Fig. [T| The As-7p transition frequency 
w 7p,4s = 12801 cm -1 corresponding to a one 781.2 nm 
photon is included in the probe pulse however, it is ex- 
cluded in the pump pulse by blocking this part of the 
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spectrum at the Fourier plane of the pulse shaper. The 
probe beam, which remains a TL pulse, is spatially over- 
lapped with the pump beam in the heated sodium cell. 

( 2) 

We evaluate the transient population P tr (i) of the 4s 
state by probing the third harmonic generation signal at 
the 7p-3s transition frequency u>7 Pl 3 S — 28541 cm^ 1 cor- 
responding to a one 259.5 nm UV-photon [29]. We as- 
sume that the probe pulse is sufficiently weak and short 
compared to the shaped pump pulse. The UV signal is 
measured within the beam propagation direction using a 
spectrometer coupled to a time-gated camera system. 

Figure [2] demonstrates experimentally the transient 
selective coherent control strategy. It presents experi- 
mental (circles) and numerical-theoretical (lines) results 
(Eqs. 3|6l for coherent control of transient two-photon 
absorption. All the data is normalized by the final pop- 
ulation P^l(t — > oo) excited by the transform- limited 
pulse. The data is presented as a function of the probe 
pulse delay. Panels (a) and (b) show different transient 
population build-ups resulting in the same TL-level final 
population. While panel (a) shows the transient pop- 
ulation for a TL pulse, i.e., both ^base^) = ^ anc ^ 

^antis = ^ ^ or an ^ w ' P ane ^ 0-0 presents the tran- 

sient population of a modified pulse based on the sym- 
metry properties of the system for which ^^^(cj) = 

remains zero for any ui, however ( w ) 1S n °t- It 

has a shape of an antisymmetric step (see Fig. [I]) with 
A = 0.55 nm and <™** sym = vr. Although the final pop- 
ulation is not affected by this phase addition the transient 
population build-up is very different as compared to the 
TL pulse. As long as for the former (a), the population 
transfer is a smooth line from zero population before the 
pulse to TL-level at the end of the pulse, for the later (b) , 
it is much longer in time and has a kink in the vicinity 
of zero probe delay. 

As another example Fig.[2](c) and (d) present the same 
scenario of selective transient coherent control for zero 
final population. Panel (c) shows the transient popula- 
tion evolution which eventually results in zero population 
transfer, i.e., so called a dark-pulse [6]. It utilizes non- 
zero base phase with a 7r-step at u^ s ^ r — 778.2 nm 
which assures zero final population. The antisymmet- 
ric addition <j> ( - ac ^ ( w ) = remains zero for any u>. 

antisym v ' J 

First, the transient population is gradually accumulated 
in the final state until it reaches a value of 20% of the 
TL-level when the probe delay is zero. Then, all the ac- 
cumulated population is transferred back to the ground 
state. In panel (c), while, the base phase with a 7r-step at 
UJ base rk ~ 772.2 nm assures a zero final population, the 

antisymmetric addition $( ac J^ ( w ) with A = 0.55 nm 
J antisym v ' 

and <^amp Vm = alters completely the transient popula- 
tion evolution without changing the zero final population. 
It has two peaks across the evolution with a deep around 
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FIG. 2: (a)-(d) Experimental (circles) and numerical- 
theoretical (lines) results for coherent control of transient two- 
photon absorption. Panels (a) and (b) show different transient 
population build-ups resulting in the same, TL-level final pop- 
ulation. Panels (c) and (d) show different transient popula- 
tion resulting in zero final population, (see text). All the 
data is normalized by P!fl(t —* oo) excited by the transform- 
limited pulse. The data is presented as a function of the probe 
pulse delay. 
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FIG. 3: Experimental (circles) and numerical-theoretical 
(lines) results for symmetry based independent control of the 
duration of the transient population evolution regardless of 
the desired final population. The left column presents the re- 
sults for TL-level final population with increasing cubic spec- 
tral phase from tpz — for the first row, ip->, = 3 ps 3 for the 
second row, ips — 5 ps 3 for the third, and if3 = 7 ps 3 for 
the forth. The same cubic phase additions were applied in 
the middle, and right rows respectively for half of the TL- 
level final population and zero final population. The data is 
presented as a function of the probe pulse delay. 



zero probe delay, and as in the case of TL-level final pop- 
ulation, in panels (a) and (b), the modified pulse induces 
much longer evolution of the transient population. It 
worth mentioning that all the experimental results are in 
excellent agrement with the theoretical calculations. 

The presented symmetry based control scheme can be 
applied to independent control of the duration of the 
transient population evolution regardless of the desired 
final population. The experimental (circles) together 
with the theoretical (line) results are presented in Fig. [3] 
In this results the final population is controlled using 
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quadratic chirp coefisient cp 3 (ps ) 

FIG. 4: Duration population evolution control. Calculated 
results show the duration of transient duration as a function 
of quadratic chirp coefficient 993 for different final populations. 



^base^) m a f° rm °f a simple 7r-step by setting the step 
at different position w^p. The duration of the transient 
population evolution is controlled using quadratic chirp 

spectral phase addition ^antfsym^ = <^3 (w - uj fg /2) 3 , 
where ip3 is a third-order coefficient. This phase addition 
is antisymmetric around the half of the two-photon tran- 
sition frequency ujf g /2. Being so, it will not influence the 
final population set by the base phase. 

The left column in Fig. [3] presents the results for TL- 
levcl final population with increasing cubic spectral phase 
from (ps = for the first row, ip^ — 3 ps 3 for the second 
row, ip 3 = 5 ps 3 for the third, and (p^ — 7 ps 3 for the 
forth. The same cubic phase additions were applied in 
the middle, and right rows respectively for half of the 
TL-level final population and zero final population. As 
Fig.[3]clearly indicates, the duration of the transient pop- 
ulation evolution is highly correlated to the magnitude of 
the quadratic chirp addition, regardless of the final pop- 
ulation. As we increase the magnitude of the quadratic 
chirp, the duration of transients also increase. This tran- 
sient elongation is very intuitive as it is well known that 
quadratic chirp produce beats in the intensity vs. time, 
causing oscillations after the main pulse. Those oscilla- 
tions become more and more pronounce as we incense the 
coefficient ^3, effectively elongating the pulse, ft worth 
mentioning that linear chirp addition is not applicable for 
the current scenario as it is symmetric around Wf g /2 so 
it will alter the final population together with the tran- 
sient duration. Although, in general all the experimen- 
tal results are in excellent agrement with the theoretical 
calculations, some discrepancies that do appear, can be 
accounted for finite pulse shaper resolution [28] and the 
high sensitivity of the signal |26j . The effect is clearly 
more pronounced for high chirp values as we get close to 
the sampling limitation of the shaper. 

To emphasize the robustness of the selective tran- 
sient coherent control which allows independent control 
of transient population evolution for a given final popula- 
tion, duration of such evolution is presented as a function 



of the quadratic chirp coefficient ip^ for different final 
populations. The calculated data for a Gaussian pulse 
corresponding to 100 fs unshaped duration is presented 
in Fig. [4] The transient duration is calculated by looking 
at 5% (relative to TL-lcvcl) deviation of the population 
from the initial and final levels. The results clearly indi- 
cate the possibility of transient duration control indepen- 
dent of the final population. The discrepancies between 
the curves for deferent final levels are accounted for the 
duration calculation procedure which does not take into 
account the relativity of the final population. The com- 
mon feature for all the curves in Fig. [4] appear to have 
little correlation for small values of the chirp and a very 
pronounced linear dependence for high chirp values. As 
for the first region, both the Q^^^u) in a form of a sim- 
ple 7r-step and the <j>( ac J^' ( w ) have similar relative ef- 
1 1 antisynr ' 

feet on the pulse duration, thus influencing the transient 
duration in the same manner. As for the linear region, 
the 7r-step influence on the pulse duration is negligible 
relative to the applied chirp, thus resulting in little im- 
pact on the duration of the transient, which is controlled 
solely by chirp magnitude. Despite its reduced correla- 
tion to the duration of the transient evolution, the basic 
part ^>b ase (w) of the applied phase is vital for control 
of the final population by setting the proper interference 
pattern. Such combination of frequency domain picture 
given in Eqs. 3]|6 with the time domain perspective gives 
a simple intuitive way to draw the control scenario for 
independent selective control of both the transient and 
final populations. While the frequency domain formula- 
tion allows us to identify the interfering pathways selec- 
tively contributing to the final and transient populations, 
the time domain picture sets the general limits on the du- 
ration of the transient evolution. 

In conclusion, full selective coherent control of tran- 
sient population is presented for the first time in which 
both transient and final populations are controlled inde- 
pendently. The control scheme is very robust, intuitive 
and general, and can be applied to any system, once its 
symmetry properties are identified. The duration of the 
transient build-up can be controlled from 100 fs to several 
picoseconds for any desired final population. Moreover, 
the presented control scheme can be extended to basically 
any structure of transient evolution by employing a semi- 
rational however automatic phase optimization based on 
closed- loop learning algorithms [30 , 3H |32l [33] . 
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